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The space-time structure of the multipion system created in central relativistic heavy- 
ion collisions is investigated. Using the microscopic transport model UrQMD we determine 
the freeze-out hypersurface from equation on pion density n(t, r) = n c . It turns out that 
for proper value of the critical energy density e c equation e(t, r) = e c gives the same freeze- 
^ . out hypersurface. It is shown that for big enough collision energies E]^ n > 40 A GeV/c 

> 8A GeV/c) the multipion system at a time moment r ceases to be one connected 
unit but splits up into two separate spatial parts (drops), which move in opposite direc- 
. tions from one another with velocities which approach the speed of light with increase 

of collision energy. This time r is approximately invariant of the collision energy, and 
P\| ' the corresponding r=const. hypersurface can serve as a benchmark for the freeze-out 

time or the transition time from the hydrostage in hybrid models. The properties of this 
hypersurface are discussed. 
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1 Introduction 

High energy heavy ion collisions at AGS, SPS, RHIC and LHC energies are well described by 
fluid dynamical models, and a substantial part of the initial beam energy is converted into 
collective flow. This collective flow is one of the most important observable phenomena arising 
from the strongly interacting, high energy density quark-gluon plasma (QGP). 

On the other hand the initial state, before local thermal and chemical equilibrium is reached 
must be described by other means. In the subsequent fluid dynamical (FD) stage in local 
equilibrium we can take advantage of an equation of state (EoS). At the final stage when the 
system becomes dilute and local thermal and chemical equilibrium cannot be maintained we 
have to use another freeze out (FO) model. Most frequently the transition from the FD to the 
FO stage is assumed to happen at a 3-dimensional space-time hypersurface where one generates 
the final out of equilibrium state by using the Cooper- Frye formalism [Tj. When applying this 
formalism one should also take care for energy, momentum and baryon charge conservation, 
as summarized recently in ref. [2]. This is also necessary if instead of immediately applying 
a post FO distribution we supplement the FD description with a non-equilibrium model, in a 
multimodule or hybrid model approach, e.g. [3]. Due to the explosive dynamics of the system 
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the FO hyper-surface has parts both with space-like and with time-like normal vectors, which 
requires a careful handling of the FO process. 

In the present work we analyze the FO hyper-surface using the ultrarelativistic quantum 
molecular dynamics (UrQMD) model [31 H], in the energy range of AGS, SPS, and RHIC 
for central heavy ion collisions and the structure of these hypersurfaces under different FO 
conditions. We assume in this work that the FO takes place in the hadronic phase and every 
hadronic species has its own FO hypersurface. Because pions represent the majority of the 
secondary particles we study first of all the FO hypersurface for the final pion gas. 



2 Calculation algorithm 

The invariant scalar n~ density is defined as [3] 

n(x) = N^(x)u^x). (1) 

Here is the particle four-flow 

N"(x) = J yLjfffap) = (n lab ,n lab vj , (2) 

where p° = + p 2 , rii^ is the density of negative pions in the laboratory frame at point 

x and v B is their average three- velocity. The quantity u^{x) is the collective four- velocity of 
negative pions. It can be written as = (7,7V), where 7(x) = 1/y/l — v 2 - The velocity v 
is tied to Local Rest frame definition. In Local Rest frame w M has only temporal component: 
u lr = (1)0). According to Eckart definition is tied to the particle (pion) flow 

u »(x) = = (7 b , 7b vJ . (3) 

(N"N V )* 

Using eq. ([1]) we can write the invariant 7r~ density according to Eckart definition as 

The invariant energy density e(x) is defined as 

e(x) = u^x) T^(x) u u {x) , (5) 
where T^ u is energy-momentum tensor 

t" u (x) = J ftjrjrffap). (6) 

To calculate quantities n(x) and e(x) we use the transport model UrQMD v2.3 (3j H]. In 
order to perform calculations we take a four-volume box with sides = 20 fm, i = t,x, y, z, 
and divide it into cells with the side length of 1 fm. For each cell the ir~ density, ni a b(x), and 
the average velocity, v B (x) of pions in that cell are calculated as a result of averaging over 1000 
collision events. Then the particle four-flow A /i (x) for this cell is calculated by using relation 
(J2J), which allows us to determine the Eckart four- velocity u^{x), see eq. Relation (pQ) then 
allows to determine the ir~ invariant density for each cell. 

In order to determine the invariant energy density in a cell, the energy-momentum tensor, 
T MV (x), is calculated in UrQMD by averaging data over 1000 events according to eq. ([6]) for 
each cell. The invariant pion energy density is then determined by eq. (jSJ). 
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UrQMD (Au+Au), E kjn =10.8A GeV 



UrQMD (Au+Au), E kin =10.8A GeV 




Figure 1: (Color online) (a) Three-dimensional hypersurface of constant invariant particle den- 
sity of negative pions for AGS conditions (E kin = 10. 8A GeV). (b) The hypersurface at constant 
invariant energy density, (c) The same hypersurface in z-r coordinates at 5.5 fm/c. 



3 Results and discussion 

In order to determine the space-time characteristics of the fireball and FO process it is use- 
ful to analyze the hypersurfaces of constant invariant ir~ particle and energy densities. Be- 
cause of the symmetry of central collisions the particle and energy densities do not depend 
on azimuthal angle if in the x-y plane when transforming to cylindrical coordinates, that is 
n(t,x,y,z) = n(t,r,z) and e(t,x,y,z) = e(t,r,z). In this case it is possible to visualize the 
constant density hypersurfaces in coordinates (t,r,z), where r = ±a/x 2 + y 2 . 

Calculation results for n(x) = n c = 0.08 /fm 3 and e(x) = e c = 0.035 GeV/fm 3 are depicted in 
Figs. dEfor central collisions with energies ranging from AGS to RHIC. It might be mentioned 
that the UrQMDv2.3 model is not known to yield accurate results at RHIC energies where 
extremely hot and dense nuclear matter is created in the early stage of collision. However, it is 
still useful to analyze the FO hypersurface within this model and to obtain at least qualitative 
space-time features of the FO process and to compare them with corresponding features at 
lower energies. 

It can be seen that the hypersurfaces, corresponding to constant particle and energy density 
virtually coincide. In fact, the value of e c = 0.035 GeV/fm 3 was specially chosen to show (to 
prove) that both FO hypersurfaces coincide with one another when they are determined by 
means of equations n(x) = n c and e(x) = e c provided by specific correspondence of n c to e c (see 
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Figs. [U15]). In case of a higher RHIC energy (^/saa = 130A GeV) the corresponding value of 
e c is higher and is equal to 0.043 GeV/fm 3 (see Fig. [6]). 




Regarding these hypersurfaces as sharp FO hypersurfaces we can determine temperature Tf 
at FO by using relation 

— = 3T f + m n — — , 7 

n c K 2 {m 7r /Ti) 

assuming that at the FO hypersurface we are still having an ideal, dilute pion gas described with 
the Boltzmann single-particle distribution, / B (p) oc exp (— \[m^ + p 2 /Tf), in the rest frame of 
the element of FO hypersurface. Here K\ and K 2 are the modified Bessel functions of the 
second kind. 

Solving the above equation yields Tf = 128 MeV for AGS, SPS and lower RHIC energies. 
For higher RHIC energy of ^saa = 130A GeV the temperature rises to 164 MeV. As in our 
calculation we assume the validity of the UrQMD dynamics both before and after FO, we have 
only one FO temperature Tf. If we would replace the pre-FO stage by fluid dynamics then the 
pre-FO temperature may be different, while the conserved currents across the front must be 
conserved. 

It is worth mentioning that the same value of temperature, calculated by the above formula, 
also appears in other space-time regions. That is, there are several hypersurfaces corresponding 
to the same value of temperature so it is only possible to define temperature at FO, but it is 
not possible to define a single, unique, connected FO hypersurface of constant temperature. 
Especially at higher energies, the constant temperature hypersurface fragments to unconnected 
pieces. 

We can define the lifetime r of the pion system (pion fireball) confined by the FO hypersur- 
face as a maximum time moment of the projection of FO hypersurface on t — z plane. It is seen 
that at AGS energies (see Fig. [Hand Tabled]) the fireball has a lifetime of about r = 9 fm/c, 
and maximum transverse and longitudinal radii of 6 and 2 fm respectively. The exact values 
of these parameters would change for different choice of n c but it can be seen that for this 
energies the longitudinal expansion of pion fireball is rather small. One qualitative feature of 
pion fireball at AGS energies is that it always stays spatially as a single freezing out object 
during the whole time evolution. 
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UrQMD (Pb+Pb), E k . =158A GeV 



UrQMD (Pb+Pb), E. =158A GeV 




Figure 3: (Color online) Same as Fig. [I] but for SPS conditions (E kin = 158A GeV). 



The situation changes at SPS energies (see Figs. [2] and Tabled]). It is seen that after some 
time td ~ 9 fm/c the pion fireball spatially breaks up into two parts (see the projection of FO 
hypersurface onto t — z plane). That is, t& is actually a fireball division time. 

Hence, at the time t = t& the fireball is separating into two pieces (drops), which lifetimes 
are counting starting from this moment. By considering the fireball just as a whole, connected 
body in the 4-dimensional space-time (i.e. a single piece of nuclear matter) we can generalize 
the definition of the fireball lifetime r. So, we define the lifetime of the fireball as that time 
span when fireball exists just as one single piece of exited nuclear matter, i.e. r = t&. After 
division of the fireball there are two drops of matter which move in opposite directions from 
one another with velocities approaching the speed of light with increase of collision energy (see 
Figs. |2H6l especially upper and lower panels of Fig. [3]). 

With the increase of energy the longitudinal radius increases while the fireball division time 
ta stays about the same. Since the velocity of pions cannot exceed the speed of light and the 
fireball formation time is very small compared to the fireball lifetime, the longitudinal radius is 
bounded from above by r. It reaches it's maximum value at SPS energies of about 7 fm, which 
is consistent with HBT interferometry measurements [6]. 

At higher energies available at RHIC (Figs. H] and Table [T]) the picture is similar to SPS 
energies: the fireball expands as a whole until about ~ 8.25 — 10 fm/c and then breaks up 
into two parts. The values of R\\ are consistent with HBT radii [8]. 

What indicates that the fireball division times, t&, and thus the fireball lifetimes, r, change 
very weakly with collision energy (see Table [1]), they vary in the range r ~ 8.25 — 10.0 fm/c. 
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UrQMD (Au+Au), s 1/2 =19.6A GeV 




20 -20 

Figure 4: (Color online) Same as Fig. 



UrQMD (Au+Au), s 1/2 =19.6A GeV 




for RHIC conditions (y/sJX = 19.64 GeV). 



Table 1: Space-time evolution parameters. 
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Figure 5: (Color online) Same as Fig. H]but for ^Jsaa — 62v4 GeV. 




The fireball space-time structure depends on the chosen value of n c (and/or, respectively 
on the value of e c ). The value of n c = 0.08 fm~ 3 was chosen so that maximum fireball size at 
SPS and RHIC energies corresponds to known interferometry data R ~ 6 — 8 fm. It can be 
argued what is the exact value of n c that corresponds to the sharp kinetic (or chemical) FO 
hypersurface; so, it is useful to investigate the structure of the fireball at different values of n c . 

One feature of the fireball structure for n c = 0.08 fm -3 is that the fireball lifetime, r 
depends weakly on collision energy (see Table [T|). The values of r for different values of n c 
and for different collisions energies are presented in Table |2j It is seen that r depends on the 
n c value (with the increase of n c the value of r decreases), but for each n c it depends very 
weakly on collision energy: for n c = 0.04 fm -3 it is in the range r ~ 12 — 14 fm/c and for 
n c = 0.12 fm -3 lifetime is in the range r ~ 6.25 — 7.5 fm/c. Thus, it can be stated that the 
fireball lifetime is approximately invariant of the collision energy. 

The hypersurfaces of constant tt~ particle and energy density can be compared to hadron 
reaction zones. Reaction zone is defined as the space-time region where a certain fraction of 
reactions of certain type took place (see Ref. [9]). The reaction zone, which contains 99% 
of all inelastic hadronic reactions (see. Fig. [7j) is related to 7r~ FO process (since inelastic 
reactions mostly involve pions) and includes the hypersurfaces of constant 7r~ particle and 
energy densities (which can be regarded as sharp FO hypersurfaces). This particular reaction 
zone is called hot fireball (see Ref. [9]). It is also seen that the reaction zone boundary (see 
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UrQMD (Pb+Pb), E kin =158A GeV 




Zone of 99% inelastic reactions 
z-r projection 
z-t projection 
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Figure 7: (Color online) The three-dimensional reaction zone, which determines the space-time 
region where 99% of all inelastic hadronic reactions for SPS conditions (-Ekm = 158A GeV) take 
place. 



Table 3: Division times of the reaction zone (different fractions of the contained 
inelastic reactions) 
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Fig. [7]) is similar to the hypersurfaces of constant particle density: at SPS and RHIC energies 
the fireball expands as a whole for some time, t, and then spatially breaks up into two parts, 
extending from each other in the ± ^-directions. 

Similarly to the case of constant particle and energy density hypersurfaces, the reaction 
zone division time depends weakly on collision energy (see Table [3]). It is seen from Table [3] 
that the values of reaction zone division time depend on the fraction of total inelastic hadronic 
reactions which are contained in that reaction zone, for instance 80% or 99% of all inelastic 
reactions. It is similar to the dependence of fireball lifetime on the value of n c (see Table [2]). Hot 
fireball (reaction zone containing 99% of inelastic hadronic reactions) division times are higher 
than corresponding values of lifetimes (division times) in Table |2] regarding constant n~ density 
hypersurfaces, therefore it corresponds to a rather low n c for negative pions. A different (lower) 
choice of percentage of inelastic reactions contained in reaction zone, for instance 80%, would 
correspond better to the constant 7r~ density hypersurfaces and their division times studied in 
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UrQMD (Au+Au), s =19.6A GeV 
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Figure 8: (Color online) The [t,z] projection of the constant 7r _ density hypersurface for (a) 
SPS (E kin = A0A GeV) and (b) RHIC (y/sJX = 19. 6 A GeV) conditions. The solid black line 
indicates space-like FO hypersurface boundary approximation according to eq. (JHJ). 



this work (see Table El division times for reaction zone containing 80% of inelastic reactions). 
From this we can conclude that the pion FO hypersurface, which is determined in accordance 
with the critical density n c = 0.08 fm~ 3 , contains the space-time volume where approximately 
80% of all inelastic hadronic reactions take place. 

Returning back to the pion FO hypersurface we reveal yet another feature: the space-like 
part (time-like normal) of the [t, z]-section of the FO hypersurface at SPS and RHIC energies 
can be approximated with a r PO = const hyperbola originating from possibly different time, 
i° , than the initial time t = 0, of the collision (see the upper boundary curve in Fig. |H]). The 
equation for this hyperbola has the following form 

*«>(*) = C + V^o + Z " ■ ( 8 ) 

The fireball lifetime, r (which is actually a fireball division time td) is the minimum time on 
the space-like FO hypersurface, i.e., r = t FO (z)\ z=Q . Therefore, if hypersurface approximation 
works well for central region, then the fireball lifetime, r, is related to hypersurface parameters 
as 

T = T + t° . 

FO po 

The [t, z] projections of constant ir~ density FO hypersurface for n c = 0.08 fm -3 for SPS 
and RHIC energies are depicted in Fig. |HJ The hypersurface boundary approximation Eq. [H] 
is depicted by solid black line. The values of parameters r FO and t° for different energies are 
presented in Table HI It is seen that the space-like boundary (time-like normal) is well approx- 
imated by hyperbola of constant proper time. Meanwhile, it turns out that the hyperbolas 
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originate from an earlier times t® Q < 0, than the initial time, t = 0, of the collision. At RHIC 
energies the originating time t® Q approaches zero as the energy rises (see Table Hj) which should 
be also the case for LHC energies. 

4 Conclusions 

The present studies based in the UrQMD calculations are able to give a reliable estimate of 
the final hadronic FO hypersurface in space-time. This is a very useful information for all 
global estimates and of the space-time development of the reaction. Furthermore it gives a 
good guidance for the selection of the FO hypersurface for multi-module or hybrid models, 
most importantly the [t, ^-section of the hypersurface. The radial section of the hypersurface 
is relatively simple to model or parametrize as shown in all figures presenting central collisions. 
We have to mention, however, that in peripheral collisions, where important flow phenomena 
are observed (e.g. the 3rd flow component or antiflow [10]) the FO hypersurface becomes 
more involved in the transverse directions also. Up to now the most frequently used direct FO 
descriptions have simplified the FO hypersurface to t = const, or r = const, hypersurfaces. 
These studies indicate that the FO hypersurface is more complex, and may require a more 
detailed work of finding such a surface as described in references [21 [Q] . We can also see that 
at high energies the time-like parts (space-like normals) of the FO hypersurface are almost 
parallel to the light cone, therefore the amount of particles freezing out through this part of 
the FO hypersurface is negligible compared to those crossing space-like part (time-like normal) 
of the FO hypersurface. In this work we also propose a simple alternative by observing that 
the space-like part (time-like normal) of the FO hypersurface can be well approximated with 
a FO-hyperbola, at r FO = const., but originating from an earlier time < 0, than the initial 
time, t = 0, of the collision. 

In several fluid dynamical models the isochronous freeze-out or isochronous transition from 
the hydro stage is assumed in hybrid models (see, e.g., [T2] ). These results can be improved 
by using an iso-r transition hypersurface [13]. The optimal choice of the proper time of the 
transition can be chosen based on the present studies (see Table 0]). As discussed above up 
to about y/sAA — 20 A GeV the origin of the freeze-out hyperbola is found to be at negative 
times, which is a consequence of the string dynamics and its hadronization in UrQMD. This 
work indicates that this option should also be taken into account at the selection of freeze-out 
or transition hypersurface in fluid dynamical and hybrid models. 

At the same time, the minimum of the FO-hyperbola r = t FO (z)\ z _ Q corresponds to the 
time moment t = r when the many-particle system ceases to be one unit but splits up into 
two separate spatial parts which move in opposite directions from one another with velocities 
approaching the speed of light with increase of collision energy. This time r is nothing more as 
a lifetime of the fireball (if we treat a fireball as one unit) and it turns out that r depends so 
weakly on collision energy that we can state that the fireball lifetime is approximately invariant 
of the collision energy. These results provide different possibilities to construct more realistic 
models for the selection of transition or freeze out hypersurfaces. 
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